Return to MUSA 801 Projects Page

This project was completed for the MUSA/Smart Cities Practicum course (MUSA 801) instructed by Michael Fichman and Matthew Harris. We are grateful to our instructors for their continued support and feedback. We would like to give special thanks to KC Filippino and Ben McFarlane from Hampton Roads Planning District Commission, and Dexter Locke from the United States Forest Service for providing data, insight, and support throughout the semester. This project would not have been possible without them.

1.Introduction

1.1 Abstract

This project aims to develop a precision forecast model for land cover change at the Chesapeake Watershed, focusing on three representative counties: Isle of Wight, James City, and Portsmouth. By leveraging high-resolution longitudinal land use and land cover data provided by the Chesapeake Conservancy, the model will predict land cover conversions from pervious to impervious surfaces. This forecast will enable land use and environmental planners to visualize and assess potential impacts on flood risk, heat stress, and heat island risks in the region. The model will be generalizable to the county level, incorporating only widely available inputs, thus allowing any municipality within the Chesapeake basin to replicate the analysis. This proof-of-concept project will demonstrate the utility of precision conservation in climate adaptation and mitigation planning and provide a valuable tool for planners and policymakers across the region.

1.2 Motivation & Use Case

To build resilient communities, the HRPDC has set up a green infrastructure plan. This plan aims to to identify and prioritize a network of valuable conservation lands in order to achieve multiple benefits, such as habitat protection, drinking water supply protection, stormwater management and recreational opportunities. A new component of this plan is to build a model for the potential future growth and identify which areas of the green infrastructure network are most at risk for development.

2. Exploratory Analysis

2.1 Understanding landcover data

The Chesapeake Conservancy supplies us with high-resolution landcover data, which is essential for our precision conservation efforts. This data is a vast raster dataset with an impressive 1-meter accuracy, offering 900 times more detail than the commonly used 30-meter resolution National Land Cover Dataset. This level of detail is critical in capturing subtle changes in land use and land cover. Within the landcover classification, pervious surfaces include categories such as tree canopy, shrub, and wetlands, which allow water to infiltrate the ground, while impervious surfaces encompass categories like roads and structures that prevent water infiltration, leading to increased runoff and potential flooding issues. This detailed classification enables us to better understand and predict land cover changes, particularly the conversion from pervious to impervious surfaces.

Data Source

2.2 Data Cleaning and Wrangling

2.3

dependent: pervious to impervious surface during the period 2014-201

independent:

lccover change 能change
1. `originallc` originallandcover

2. `lcp` permeable/impermeable:
< 6 - 0
>= 6 - 1

3. `lcchange` whether landcover has changed: 1

4. `lc` 3*3的permeable rate:
0-1(0.1,0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8)

5. `popchange` per m^2

6. `pctwhitechange`

7. `unitchange` per m^2

8. `medhhincchange` per m^2

9. `road` road 1/0

10. `water` water

11. `canopy` tree canopy:
tree canopy - 1
tree canopy over structure -0.5
tree canopy over other -0.5
tree canopy over road - 0.5
other - 0

12. `perm` permeable:
2 -shrub - 0.3
4 - herb - 0.3
5 - wetand - 0.4
other -0

13. `barren` barren:
6 - barren - 1
not barren - 0

14. `struct` structure:
7 - structure - 1
not structure - 0

15. `other` other:
8 - other - 1
not other - 0

16. `terrain`dem

17. `slope`slope

18. `area` m^2 per block group

19. `geoid`-- blockgroup id

20. `soil type`-

模型encode:categorial:尤其是

possible features per 10*10 :

Road density distance to rivers

impervious / pervious landcover type population density population change per 4 year median hh income median hh income change per 4 year percentage of white population percentage of white population change per 4 year percentage of changing from impervious to pervious/from pervious to impervious in 4 years predicted landcover type:

possible features per block group :

percentage of impervious / pervious change: population population change per 4 year median hh income median hh income change per 4 year percentage of white population percentage of white population change per 4 year percentage of changing from impervious to pervious/from pervious to impervious in 4 years predicted landcover type composition:

possible feautres per county:

percentage of impervious / pervious change: population population change per 4 year median hh income median hh income change per 4 year percentage of white population percentage of white population change per 4 year percentage of changing from impervious to pervious/from pervious to impervious in 4 years predicted landcover type composition: development speed defined by: speed of population change and landcover change speed Agricultural land use, vegetation cover and water loss:

LS0tDQp0aXRsZTogIlByZWNpc2lvbiBGb3JlY2FzdHMgb2YgTGFuZCBDb3ZlciBDaGFuZ2UiDQpzdWJ0aXRsZTogIkNoZXNhcGVha2UgV2F0ZXJzaGVkIg0KYXV0aG9yOiAiWXVld2VuIERhaSwgU2h1amluZyBZaSwgWGluZ2UgWmhhbmciDQpkYXRlOiAiMjAyMy0wNC0yNSINCm91dHB1dDogDQogIGh0bWxfZG9jdW1lbnQ6DQogICAgdG9jOiB0cnVlDQogICAgdG9jX2Zsb2F0OiB0cnVlDQogICAgY29kZV9mb2xkaW5nOiAiaGlkZSINCiAgICBjb2RlX2Rvd25sb2FkOiB0cnVlDQotLS0NCg0KYGBge3Igc2V0dXAsIGluY2x1ZGU9RkFMU0V9DQprbml0cjo6b3B0c19jaHVuayRzZXQoZWNobyA9IFRSVUUpDQpgYGANCg0KW1JldHVybiB0byBNVVNBIDgwMSBQcm9qZWN0cyBQYWdlXShodHRwczovL3Blbm5tdXNhLmdpdGh1Yi5pby9NVVNBXzgwMS5pby8pDQoNClRoaXMgcHJvamVjdCB3YXMgY29tcGxldGVkIGZvciB0aGUgTVVTQS9TbWFydCBDaXRpZXMgUHJhY3RpY3VtIGNvdXJzZSAoTVVTQSA4MDEpIGluc3RydWN0ZWQgYnkgTWljaGFlbCBGaWNobWFuIGFuZCBNYXR0aGV3IEhhcnJpcy4gV2UgYXJlIGdyYXRlZnVsIHRvIG91ciBpbnN0cnVjdG9ycyBmb3IgdGhlaXIgY29udGludWVkIHN1cHBvcnQgYW5kIGZlZWRiYWNrLiBXZSB3b3VsZCBsaWtlIHRvIGdpdmUgc3BlY2lhbCB0aGFua3MgdG8gS0MgRmlsaXBwaW5vIGFuZCBCZW4gTWNGYXJsYW5lIGZyb20gSGFtcHRvbiBSb2FkcyBQbGFubmluZyBEaXN0cmljdCBDb21taXNzaW9uLCBhbmQgRGV4dGVyIExvY2tlIGZyb20gdGhlIFVuaXRlZCBTdGF0ZXMgRm9yZXN0IFNlcnZpY2UgZm9yIHByb3ZpZGluZyBkYXRhLCBpbnNpZ2h0LCBhbmQgc3VwcG9ydCB0aHJvdWdob3V0IHRoZSBzZW1lc3Rlci4gVGhpcyBwcm9qZWN0IHdvdWxkIG5vdCBoYXZlIGJlZW4gcG9zc2libGUgd2l0aG91dCB0aGVtLg0KDQojIyAxLkludHJvZHVjdGlvbg0KIyMjIDEuMSBBYnN0cmFjdA0KVGhpcyBwcm9qZWN0IGFpbXMgdG8gZGV2ZWxvcCBhIHByZWNpc2lvbiBmb3JlY2FzdCBtb2RlbCBmb3IgbGFuZCBjb3ZlciBjaGFuZ2UgYXQgdGhlIENoZXNhcGVha2UgV2F0ZXJzaGVkLCBmb2N1c2luZyBvbiB0aHJlZSByZXByZXNlbnRhdGl2ZSBjb3VudGllczogSXNsZSBvZiBXaWdodCwgSmFtZXMgQ2l0eSwgYW5kIFBvcnRzbW91dGguIEJ5IGxldmVyYWdpbmcgaGlnaC1yZXNvbHV0aW9uIGxvbmdpdHVkaW5hbCBsYW5kIHVzZSBhbmQgbGFuZCBjb3ZlciBkYXRhIHByb3ZpZGVkIGJ5IHRoZSBDaGVzYXBlYWtlIENvbnNlcnZhbmN5LCB0aGUgbW9kZWwgd2lsbCBwcmVkaWN0IGxhbmQgY292ZXIgY29udmVyc2lvbnMgZnJvbSBwZXJ2aW91cyB0byBpbXBlcnZpb3VzIHN1cmZhY2VzLiBUaGlzIGZvcmVjYXN0IHdpbGwgZW5hYmxlIGxhbmQgdXNlIGFuZCBlbnZpcm9ubWVudGFsIHBsYW5uZXJzIHRvIHZpc3VhbGl6ZSBhbmQgYXNzZXNzIHBvdGVudGlhbCBpbXBhY3RzIG9uIGZsb29kIHJpc2ssIGhlYXQgc3RyZXNzLCBhbmQgaGVhdCBpc2xhbmQgcmlza3MgaW4gdGhlIHJlZ2lvbi4gVGhlIG1vZGVsIHdpbGwgYmUgZ2VuZXJhbGl6YWJsZSB0byB0aGUgY291bnR5IGxldmVsLCBpbmNvcnBvcmF0aW5nIG9ubHkgd2lkZWx5IGF2YWlsYWJsZSBpbnB1dHMsIHRodXMgYWxsb3dpbmcgYW55IG11bmljaXBhbGl0eSB3aXRoaW4gdGhlIENoZXNhcGVha2UgYmFzaW4gdG8gcmVwbGljYXRlIHRoZSBhbmFseXNpcy4gVGhpcyBwcm9vZi1vZi1jb25jZXB0IHByb2plY3Qgd2lsbCBkZW1vbnN0cmF0ZSB0aGUgdXRpbGl0eSBvZiBwcmVjaXNpb24gY29uc2VydmF0aW9uIGluIGNsaW1hdGUgYWRhcHRhdGlvbiBhbmQgbWl0aWdhdGlvbiBwbGFubmluZyBhbmQgcHJvdmlkZSBhIHZhbHVhYmxlIHRvb2wgZm9yIHBsYW5uZXJzIGFuZCBwb2xpY3ltYWtlcnMgYWNyb3NzIHRoZSByZWdpb24uDQoNCg0KIVtdKGltYWdlcy9IUlBEQy5qcGcpDQoNCg0KIyMjIDEuMiBNb3RpdmF0aW9uICYgVXNlIENhc2UNCg0KVG8gYnVpbGQgcmVzaWxpZW50IGNvbW11bml0aWVzLCB0aGUgSFJQREMgaGFzIHNldCB1cCBhIGdyZWVuIGluZnJhc3RydWN0dXJlIHBsYW4uIFRoaXMgcGxhbiBhaW1zIHRvIHRvIGlkZW50aWZ5IGFuZCBwcmlvcml0aXplIGEgbmV0d29yayBvZiB2YWx1YWJsZSBjb25zZXJ2YXRpb24gbGFuZHMgaW4gb3JkZXIgdG8gYWNoaWV2ZSBtdWx0aXBsZSBiZW5lZml0cywgc3VjaCBhcyBoYWJpdGF0IHByb3RlY3Rpb24sIGRyaW5raW5nIHdhdGVyIHN1cHBseSBwcm90ZWN0aW9uLCBzdG9ybXdhdGVyIG1hbmFnZW1lbnQgYW5kIHJlY3JlYXRpb25hbCBvcHBvcnR1bml0aWVzLiBBIG5ldyBjb21wb25lbnQgb2YgdGhpcyBwbGFuIGlzIHRvIGJ1aWxkIGEgbW9kZWwgZm9yIHRoZSBwb3RlbnRpYWwgZnV0dXJlIGdyb3d0aCBhbmQgaWRlbnRpZnkgd2hpY2ggYXJlYXMgb2YgdGhlIGdyZWVuIGluZnJhc3RydWN0dXJlIG5ldHdvcmsgYXJlIG1vc3QgYXQgcmlzayBmb3IgZGV2ZWxvcG1lbnQuDQoNCiMjIDIuIEV4cGxvcmF0b3J5IEFuYWx5c2lzDQojIyMgMi4xIFVuZGVyc3RhbmRpbmcgbGFuZGNvdmVyIGRhdGENClRoZSBDaGVzYXBlYWtlIENvbnNlcnZhbmN5IHN1cHBsaWVzIHVzIHdpdGggaGlnaC1yZXNvbHV0aW9uIGxhbmRjb3ZlciBkYXRhLCB3aGljaCBpcyBlc3NlbnRpYWwgZm9yIG91ciBwcmVjaXNpb24gY29uc2VydmF0aW9uIGVmZm9ydHMuIFRoaXMgZGF0YSBpcyBhIHZhc3QgcmFzdGVyIGRhdGFzZXQgd2l0aCBhbiBpbXByZXNzaXZlIDEtbWV0ZXIgYWNjdXJhY3ksIG9mZmVyaW5nIDkwMCB0aW1lcyBtb3JlIGRldGFpbCB0aGFuIHRoZSBjb21tb25seSB1c2VkIDMwLW1ldGVyIHJlc29sdXRpb24gTmF0aW9uYWwgTGFuZCBDb3ZlciBEYXRhc2V0LiBUaGlzIGxldmVsIG9mIGRldGFpbCBpcyBjcml0aWNhbCBpbiBjYXB0dXJpbmcgc3VidGxlIGNoYW5nZXMgaW4gbGFuZCB1c2UgYW5kIGxhbmQgY292ZXIuIFdpdGhpbiB0aGUgbGFuZGNvdmVyIGNsYXNzaWZpY2F0aW9uLCBwZXJ2aW91cyBzdXJmYWNlcyBpbmNsdWRlIGNhdGVnb3JpZXMgc3VjaCBhcyB0cmVlIGNhbm9weSwgc2hydWIsIGFuZCB3ZXRsYW5kcywgd2hpY2ggYWxsb3cgd2F0ZXIgdG8gaW5maWx0cmF0ZSB0aGUgZ3JvdW5kLCB3aGlsZSBpbXBlcnZpb3VzIHN1cmZhY2VzIGVuY29tcGFzcyBjYXRlZ29yaWVzIGxpa2Ugcm9hZHMgYW5kIHN0cnVjdHVyZXMgdGhhdCBwcmV2ZW50IHdhdGVyIGluZmlsdHJhdGlvbiwgbGVhZGluZyB0byBpbmNyZWFzZWQgcnVub2ZmIGFuZCBwb3RlbnRpYWwgZmxvb2RpbmcgaXNzdWVzLiBUaGlzIGRldGFpbGVkIGNsYXNzaWZpY2F0aW9uIGVuYWJsZXMgdXMgdG8gYmV0dGVyIHVuZGVyc3RhbmQgYW5kIHByZWRpY3QgbGFuZCBjb3ZlciBjaGFuZ2VzLCBwYXJ0aWN1bGFybHkgdGhlIGNvbnZlcnNpb24gZnJvbSBwZXJ2aW91cyB0byBpbXBlcnZpb3VzIHN1cmZhY2VzLg0KDQohW10oaW1hZ2VzL0RhdGEucG5nKQ0KDQpbRGF0YSBTb3VyY2VdKGh0dHBzOi8vd3d3LmNoZXNhcGVha2Vjb25zZXJ2YW5jeS5vcmcvY29uc2VydmF0aW9uLWlubm92YXRpb24tY2VudGVyL2hpZ2gtcmVzb2x1dGlvbi1kYXRhL2x1bGMtZGF0YS1wcm9qZWN0LTIwMjIvKQ0KDQoNCiMjIyAyLjIgRGF0YSBDbGVhbmluZyBhbmQgV3JhbmdsaW5nDQoNCiMjIyAyLjMgDQoNCiMjIyMgZGVwZW5kZW50OiBwZXJ2aW91cyB0byBpbXBlcnZpb3VzIHN1cmZhY2UgZHVyaW5nIHRoZSBwZXJpb2QgMjAxNC0yMDENCiMjIyMgaW5kZXBlbmRlbnQ6DQogICAgbGNjb3ZlciBjaGFuZ2Ug6IO9Y2hhbmdlDQogICAgMS4gYG9yaWdpbmFsbGNgIG9yaWdpbmFsbGFuZGNvdmVyDQogICAgDQogICAgMi4gYGxjcGAgcGVybWVhYmxlL2ltcGVybWVhYmxlOg0KICAgIDwgNiAtIDANCiAgICA+PSA2IC0gMQ0KICAgIA0KICAgIDMuIGBsY2NoYW5nZWAgd2hldGhlciBsYW5kY292ZXIgaGFzIGNoYW5nZWQ6IDENCiAgICANCiAgICA0LiBgbGNgIDMqM+eahHBlcm1lYWJsZSByYXRlOg0KICAgIDAtMe+8iDAuMe+8jDAuMu+8jCAwLjPvvIwgMC4077yMIDAuNe+8jCAwLjbvvIwgMC4377yMIDAuOO+8iQ0KICAgIA0KICAgIDUuIGBwb3BjaGFuZ2VgIHBlciBtXjINCiAgICANCiAgICA2LiBgcGN0d2hpdGVjaGFuZ2VgDQogICAgDQogICAgNy4gYHVuaXRjaGFuZ2VgIHBlciBtXjINCiAgICANCiAgICA4LiBgbWVkaGhpbmNjaGFuZ2VgIHBlciBtXjINCiAgICANCiAgICA5LiBgcm9hZGAgcm9hZCAxLzANCiAgICANCiAgICAxMC4gYHdhdGVyYCB3YXRlcg0KICAgIA0KICAgIDExLiBgY2Fub3B5YCB0cmVlIGNhbm9wee+8mg0KICAgIHRyZWUgY2Fub3B5IC0gMQ0KICAgIHRyZWUgY2Fub3B5IG92ZXIgc3RydWN0dXJlIC0wLjUNCiAgICB0cmVlIGNhbm9weSBvdmVyIG90aGVyIC0wLjUNCiAgICB0cmVlIGNhbm9weSBvdmVyIHJvYWQgLSAwLjUNCiAgICBvdGhlciAtIDANCiAgICANCiAgICAxMi4gYHBlcm1gIHBlcm1lYWJsZToNCiAgICAyIC1zaHJ1YiAtIDAuMw0KICAgIDQgLSBoZXJiIC0gMC4zDQogICAgNSAtIHdldGFuZCAtIDAuNA0KICAgIG90aGVyIC0wDQogICAgDQogICAgMTMuIGBiYXJyZW5gIGJhcnJlbjoNCiAgICA2IC0gYmFycmVuIC0gMQ0KICAgIG5vdCBiYXJyZW4gLSAwDQogICAgDQogICAgMTQuIGBzdHJ1Y3RgIHN0cnVjdHVyZToNCiAgICA3IC0gc3RydWN0dXJlIC0gMQ0KICAgIG5vdCBzdHJ1Y3R1cmUgLSAwDQogICAgDQogICAgMTUuIGBvdGhlcmAgb3RoZXI6DQogICAgOCAtIG90aGVyIC0gMQ0KICAgIG5vdCBvdGhlciAtIDANCiAgICANCiAgICAxNi4gYHRlcnJhaW5gZGVtDQogICAgDQogICAgMTcuIGBzbG9wZWBzbG9wZQ0KICAgIA0KICAgIDE4LiBgYXJlYWAgbV4yIHBlciBibG9jayBncm91cA0KICAgIA0KICAgIDE5LiBgZ2VvaWRgLS0gYmxvY2tncm91cCBpZA0KICAgIA0KICAgIDIwLiBgc29pbCB0eXBlYC0NCg0KDQrmqKHlnotlbmNvZGXvvJpjYXRlZ29yaWFs77ya5bCk5YW25pivDQoNCnBvc3NpYmxlIGZlYXR1cmVzIHBlciAxMCoxMCA6DQoNClJvYWQgZGVuc2l0eQ0KKipkaXN0YW5jZSB0byByaXZlcnMqKg0KDQppbXBlcnZpb3VzIC8gcGVydmlvdXMNCmxhbmRjb3ZlciB0eXBlDQpwb3B1bGF0aW9uIGRlbnNpdHkNCnBvcHVsYXRpb24gY2hhbmdlIHBlciA0IHllYXINCm1lZGlhbiBoaCBpbmNvbWUNCm1lZGlhbiBoaCBpbmNvbWUgY2hhbmdlIHBlciA0IHllYXINCnBlcmNlbnRhZ2Ugb2Ygd2hpdGUgcG9wdWxhdGlvbg0KcGVyY2VudGFnZSBvZiB3aGl0ZSBwb3B1bGF0aW9uIGNoYW5nZSBwZXIgNCB5ZWFyDQpwZXJjZW50YWdlIG9mIGNoYW5naW5nIGZyb20gaW1wZXJ2aW91cyB0byBwZXJ2aW91cy9mcm9tIHBlcnZpb3VzIHRvIGltcGVydmlvdXMgaW4gNCB5ZWFycw0KcHJlZGljdGVkIGxhbmRjb3ZlciB0eXBlOg0KDQpwb3NzaWJsZSBmZWF0dXJlcyBwZXIgYmxvY2sgZ3JvdXAgOg0KDQpwZXJjZW50YWdlIG9mIGltcGVydmlvdXMgLyBwZXJ2aW91cyBjaGFuZ2U6DQpwb3B1bGF0aW9uDQpwb3B1bGF0aW9uIGNoYW5nZSBwZXIgNCB5ZWFyDQptZWRpYW4gaGggaW5jb21lDQptZWRpYW4gaGggaW5jb21lIGNoYW5nZSBwZXIgNCB5ZWFyDQpwZXJjZW50YWdlIG9mIHdoaXRlIHBvcHVsYXRpb24NCnBlcmNlbnRhZ2Ugb2Ygd2hpdGUgcG9wdWxhdGlvbiBjaGFuZ2UgcGVyIDQgeWVhcg0KcGVyY2VudGFnZSBvZiBjaGFuZ2luZyBmcm9tIGltcGVydmlvdXMgdG8gcGVydmlvdXMvZnJvbSBwZXJ2aW91cyB0byBpbXBlcnZpb3VzIGluIDQgeWVhcnMNCnByZWRpY3RlZCBsYW5kY292ZXIgdHlwZSBjb21wb3NpdGlvbjoNCg0KcG9zc2libGUgZmVhdXRyZXMgcGVyIGNvdW50eToNCg0KcGVyY2VudGFnZSBvZiBpbXBlcnZpb3VzIC8gcGVydmlvdXMgY2hhbmdlOg0KcG9wdWxhdGlvbg0KcG9wdWxhdGlvbiBjaGFuZ2UgcGVyIDQgeWVhcg0KbWVkaWFuIGhoIGluY29tZQ0KbWVkaWFuIGhoIGluY29tZSBjaGFuZ2UgcGVyIDQgeWVhcg0KcGVyY2VudGFnZSBvZiB3aGl0ZSBwb3B1bGF0aW9uDQpwZXJjZW50YWdlIG9mIHdoaXRlIHBvcHVsYXRpb24gY2hhbmdlIHBlciA0IHllYXINCnBlcmNlbnRhZ2Ugb2YgY2hhbmdpbmcgZnJvbSBpbXBlcnZpb3VzIHRvIHBlcnZpb3VzL2Zyb20gcGVydmlvdXMgdG8gaW1wZXJ2aW91cyBpbiA0IHllYXJzDQpwcmVkaWN0ZWQgbGFuZGNvdmVyIHR5cGUgY29tcG9zaXRpb246DQpkZXZlbG9wbWVudCBzcGVlZCBkZWZpbmVkIGJ5OiBzcGVlZCBvZiBwb3B1bGF0aW9uIGNoYW5nZSBhbmQgbGFuZGNvdmVyIGNoYW5nZSBzcGVlZA0KQWdyaWN1bHR1cmFsIGxhbmQgdXNlLCB2ZWdldGF0aW9uIGNvdmVyIGFuZCB3YXRlciBsb3NzOg0KDQoNCg0K